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\Q , An international search is currently underway for the quark-gluon plasma — a pre- 

dicted new phase of nuclear matter where quarks roam almost freely throughout the 
medium instead of being confined to individual nucleons. 1 ' 2 Such a plasma could be 
formed through the compression and excitation that occur when nuclei collide at rel- 
ativistic speeds. With increasing compression the nucleons overlap sufficiently that 
they should lose their individual identity and transform into deconfined quarks, and 
with increasing excitation the many pions that are produced overlap sufficiently that 
they should lose their individual identity and transform into deconfined quarks and 



ant i- quarks 



Experimental identification of the quark-gluon plasma, as well as understanding 
other aspects of the process, will require knowing the overall spacetime evolution of 
the hot, dense hadronic matter that is produced in relativistic heavy- ion collisions. 
The spacetime evolution of this hadronic matter can in principle be extracted from 
experimental measurements of invariant one-particle multiplicity distributions and two- 
particle correlations in emitted pions, kaons, and other particles. The foundations 
for two-particle correlations were laid in the 1950s by Hanbury Brown and Twiss, 3 
who used two-photon correlations to measure the size of stars, and by Goldhaber et 
al., 4 who used two-pion correlations to measure the size of the interaction region in 
antiproton annihilation. Following this pioneering work, many researchers have already 
analyzed correlations among pions and among kaons produced in relativistic heavy-ion 
collisions in terms of simple models to obtain some limited information about the size 
and duration of the emitting source. However, because of the simplicity and/or lack of 
covariance of the models that have been used, the spatial and time extensions of the 
emitting source resulting from these analyses have frequently been intertangled, and 
most of the presently available results may therefore be regarded as exploratory. 
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SOURCE MODEL 



We introduce here a new realistic expanding source model for invariant one-particle 
multiplicity distributions and two-particle correlations in nearly central relativistic 
heavy-ion collisions that contains nine adjustable parameters, which are necessary and 
sufficient to properly characterize the gross properties of the source during its freezeout 
from a hydrodynamical fluid into a collection of noninteracting, free-streaming hadrons. 
These nine physically relevant parameters fall into three categories of three parame- 
ters each, with the first category corresponding to the source's longitudinal motion, 
the second category corresponding to its transverse motion, and the third category 
corresponding to its intrinsic properties. 

The three longitudinal parameters are the rapidity y s of the source's center relative 
to the laboratory frame (in terms of which the velocity v s of the source's center relative 
to the laboratory frame is given by v s = tanh?/ s ), the longitudinal spacetime rapidity 
i] of the right-hand end of the source in its own frame (in terms of which the velocity 
V£ of the right-hand end of the source in its own frame is given by vi = tanhr^o), and 
the longitudinal freezeout proper time T{ (in terms of which the longitudinal radius at 
the end of freezeout is given by Re = Tf sinh^o)- We assume that the source is boost 
invariant within the limited region between its two ends, 5,6 and that it starts expanding 
from an infinitesimally thin disk at time t = 0. 

The three transverse parameters are the transverse velocity v t and transverse radius 
-Rt of the source at the beginning of freezeout and a transverse freezeout coefficient at 
that is related to the width At in proper time during which freezeout occurs and that 
determines the shape of the freezeout hypersurface. The transverse velocity at any point 
on the freezeout hypersurface is assumed to be linear in the transverse coordinate p. 
As illustrated in Fig. [I] for the reaction considered here, 7-9 freezeout proceeds inward 
from the initial point p — Rt, z — to the source's center and then to the source's ends. 



Figure 1. Freezeout hypersurface, which specifies the positions in spacetime where the expanding 
hydrodynamical fluid is converted into a collection of noninteracting, free-streaming hadrons. 
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In our model, the coordinate-space freezeout surface is a hyperboloid of revolution of 
first one sheet and then of two sheets, which moves through the expanding source 
like the collapsing neck of a fissioning nucleus in the three-quadratic-surface shape 
parameterization. 10 

The three intrinsic parameters are the nuclear temperature T, the ratio fib/T of 
the baryon chemical potential to the temperature (from which /ib itself can be readily 
calculated), and the fraction A,,- of pions that are produced incoherently. (The analogous 
quantity for kaons is held fixed at unity, in accordance with theoretical expectations.) 

For a particular type of particle, the invariant one-particle multiplicity distribution 
and two-particle correlation are calculated in terms of a Wigner distribution function 
that includes both a direct term 11 and a term corresponding to 10 resonance decays, 12 
namely the decay of meson resonances with masses below 900 MeV and strongly de- 
caying baryon resonances with masses below 1410 MeV. Integration of the direct part 
of the Wigner distribution function over spacetime leads to the Cooper-Frye formula 
for the invariant one-particle multiplicity distribution. 13 



APPLICATION TO NEARLY CENTRAL Si + Au COLLISIONS 

As an initial application, we apply our model to the analysis of invariant 7r + , 
7r~, K + , and K~ one-particle multiplicity distributions 7,8 and 7r + and K + two-particle 
correlations 9 for nearly central Si + Au collisions at pi a b/^4 = 14.6 GeV/c. Figure [2] 



Figure 2. Comparison between model predictions and experimental data ' for the invariant 7r 
one-particle multiplicity distribution. Results for successively increasing values of the particle 
rapidity y p relative to the laboratory frame are multiplied by 10 _1 for visual separation. 
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Figure 3. Dependence of the predicted K + two-particle correlation function C upon the 
longitudinal and "out" momentum differences, for fixed values of the other three quantities upon 
which C depends. 

shows a comparison between model predictions and experimental data for an invariant 
one-particle multiplicity distribution, and Fig. |] shows the dependence of a two-particle 
correlation function upon the longitudinal and "out" momentum differences. 14 

The nine adjustable parameters of our model are determined by minimizing x 2 
with a total of 1416 data points for the six types of data considered, so the number 
of degrees of freedom v is 1407. The error for each point is calculated as the square 
root of the sum of the squares of its statistical error and its systematic error, with a 
systematic error of 15% for 7r + , 7r _ , and K + one-particle multiplicity distributions, 20% 
for the K~ one-particle multiplicity distribution, and zero for tc + and K + two-particle 
correlations. 7-9 The resulting value of x 2 is 1484.6, which corresponds to an acceptable 
value of x 2 jv = 1.055. The values of the nine parameters determined this way, along 
with their uncertainties at 99% confidence limits on all nine parameters considered 
jointly, are given in Table 1. 

Table 1. Nine adjusted source freezeout parameters. 

Value and uncertainty 
Parameter at 99% confidence 

Source rapidity y s 
Longitudinal spacetime rapidity 770 
Longitudinal freezeout proper time Tf 
Transverse freezeout velocity Vt 
Transverse freezeout radius Rt 
Transverse freezeout coefficient a t 
Nuclear temperature T 
Baryon chemical potential fJ-h/T 
Pion incoherence fraction A w 



1.355 ± 0.066 

1.47 ± 0.13 
8.2 ± 2.2 fm/c 
0.683 ± 0.048 c 
8.0 ± 1.7 fm 
-0.86 t H 
92.9 ± 4.4 MeV 
5.97 ± 0.56 
0.65 ± 0.11 
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Table 2. Additional calculated freezeout quantities. 



Quantity 



Value and uncertainty 
at 99% confidence 



Source velocity v s 
Longitudinal velocity vg 
Longitudinal freezeout radius Re 
Beginning freezeout time t\ 
Freezeout time ti at source center 
Final freezeout time £3 
Freezeout width At in proper time a 
Baryon chemical potential /Zb 
Strangeness chemical potential fi s 
Isospin chemical potential /ij 

Number -B pro j of baryons originating from projectile 
Number B tar of baryons originating from target 
Total number B to t of baryons in source 
Baryon density n\ at beginning of freezeout 13 
Baryon density n s along symmetry axis 



16.9 t 5 4 6 g fm 
3.1 tH fm/c 
8.2 ± 2.2 fm/c 
18.8 tH fm/ c 
5.9 ttt fm /c 
554 tH MeV 
75 t\l MeV 
-5.3 MeV 
26.1 tU 




0.057 tg° 032 fm" 3 
0222 + 0097 fm" 3 

u - uzzz -0.0069 1111 




Calculated under the additional assumption that the exterior matter at z = 
that freezes out first has been moving with constant transverse velocity v t 
from time t = until time t±. 

'The upper limit of 00 for this quantity arises because the beginning freezeout 
time ti could be zero, at which time the shape is an infinitesimally thin disk. 



From these underlying nine adjustable parameters we are able to calculate several 
additional freezeout quantities of physical interest and their uncertainties at 99% con- 
fidence limits, as shown in Table 2. For this reaction, the source-frame longitudinal 
velocity V£ is 0.900 to!o29 c an d the transverse freezeout velocity v t is 0.683 ± 0.048 c. 
These relatively large values mean that the contribution to the energy of the produced 
pions and kaons from collective flow is substantial, which leads to a moderately low 
freezeout nuclear temperature T of 92.9 ± 4.4 MeV. The transverse radius R t at the 
beginning of freezeout is 8.0 ±1.7 fm, and the longitudinal radius at the end of 
freezeout is 16.9 t^J fm. 

We also investigated the dependence of x 2 upon nuclear temperature T in the 
interval 50 to 250 MeV, with the remaining eight parameters varied. On both sides 
of the absolute minimum at T = 92.9 MeV, x 2 rises steeply. However, for values 
of T greater than 129 MeV, the minimum solution switches to an unphysical branch 
corresponding to the nearly instantaneous freezeout of a thin disk of large radius. 

The predictive power of the model can be ascertained by comparing calculated and 
experimental quantities that were not included in the adjustment. For example, the 
number -B pro j of baryons in the source that originate from the projectile is calculated 
to be 26.1 tjjj, which agrees well with the value 28 corresponding to a Si projectile. 
In addition, although proton data were not used in our adjustment, the calculated 
invariant proton one-particle multiplicity distribution agrees moderately well with the 
experimental distribution for particle rapidities relative to the laboratory frame of 1.3 
or greater. For lower rapidities, however, the data show far more low-p t protons than 
are predicted. We regard these excess protons as target spectators that have interacted 
somewhat, but not enough to be considered part of the expanding system. 
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CONCLUSIONS 



We have introduced a new realistic expanding source model containing the neces- 
sary and sufficient number of adjustable parameters for properly describing invariant 
one-particle multiplicity distributions and two-particle correlations in nearly central 
relativistic heavy-ion collisions. On the basis of this model, we have extracted the 
freezeout properties of the the hot, dense hadronic matter that is produced in nearly 
central Si + Au collisions at pi a bA4 = 14.6 GeV by minimizing \ 2 with a total of 1416 
data points for invariant tt + , tt~ , K + , and K~ one-particle multiplicity distributions 
and 7r + and K + two-particle correlations. The extracted properties include not only the 
nine adjustable parameters contained in the model but also several additional freezeout 
quantities that can be calculated in terms of the underlying nine parameters. In all 
cases, we have presented uncertainties at 99% confidence limits. 

For the reaction considered here, the extracted freezeout nuclear temperature T 
is 92.9 ± 4.4 MeV, the local-rest-frame baryon density n s along the symmetry axis is 
0.0222 lo!oo69 f m ~ 3 ; an d the longitudinal freezeout proper time T{ is 8.2 ± 2.2 fm/c. 
We hope that our model will be used in the future to systematically analyze other 
existing and forthcoming data on invariant one-particle multiplicity distributions and 
two-particle correlations to extract freezeout properties as functions of bombarding 
energy and target-projectile combinations. A sharp discontinuity in the dependence of 
the extracted freezeout properties upon these quantities could signal the onset of the 
long-sought-for quark-gluon plasma or other new physics. 
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Figure 2 
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Figure 3 



